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""■■■^ ■ Abstract 
43 ' 

We investigate numerically models of the static spherically symmetric boson-fermion stars in 
the scalar-tensor theory of gravity with massive dilaton field. The proper mathematical model of 
such stars is interpreted as a nonlinear two-parametric eigenvalue problem with unknown internal 
boundary. To solve this problem the Continuous Analogue of Newton Method is used. 

Keywords: mixed fermion-boson stars, scalar-tensor theory of gravity, massive dilaton field, 
two-parametric nonlinear eigenvalue problem, Continuous Analog of the Newton Method. 
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H : 1 Main model 

cz i 

Boson stars are gravitationally bound macroscopic quantum states made up of scalar bosons. They 
differ from the usual fermionic stars in that they are only prevented from collapsing gravitationally 
by the Heisenberg uncertainty principle. For a self-interacting boson field the mass of the boson star, 
even for small values of the coupling constant, turns out to be of the order of Chandrasekhar's mass 
when the boson mass is similar to proton mass. Thus, the boson stars arise as possible candidates 
for non-baryonic dark matter in the universe and, consequently, as a possible solution of one of the 
outstanding problems in modern astrophysics: the problem of non luminous matter in the universe. 
Most of the stars are of primordial origin, being formed from an original gas of fermions and bosons in 
the early universe. That is why it should be expected that most stars are a mixture of both fermions 
and bosons in different proportions. 

Boson-fermion stars are also a good model to understand more about the nature of strong gravita- 
tional fields not only in general relativity but also in the other theories of gravity. 

The most natural and promising generalizations of general relativity are the scalar-tensor theories 
of gravity || . In these theories the gravity is mediated not only by a tensor field (the metric of space- 
time) but also by a scalar field (the dilaton). These dilatonic theories of gravity contain arbitrary 
functions of the scalar field that determine the gravitational "constant" as a dynamical variable and 
the strength of the coupling between the scalar field and matter. It should be stressed that specific 

"Talk given by S.S. Yazadjiev and T.L. Boyadjiev at the 4 th General Conference of the Balkan Physical Union, Veliko 
Turnovo , Bulgaria, 22-25 August, 2000 



1 



scalar-tensor theories of gravity arise naturally as a low energy limit of the string theory jj] which is 
the most promising modern model of unification of all fundamental physical interactions. 

Boson stars in the scalar-tensor theories of gravity with massless dilaton have been widely inves- 
tigated recently both numerically and analytically || (and references therein). Mixed boson- fermion 
stars in scalar tensor theories of gravity however have not been investigated so far in contrast to general 
relativistic case where boson- fermion stars have been investigated ||. 

We will consider a static and spherically symmetric mixed boson-fermion star in asymptotic flat 
space-time. Then the main system of differential equations can be written in the following form |j| : 
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Here the metric functions A(r) and u(r), the dilaton field ip(r), the densities of the bosonic a(r) 
and the fermionic /i(r) matter are unknown functions of the radial coordinate r € [0,i? s ] U (R s ,oo), 
where R s is the unknown radius of the star. The quantity Q (the time-frequency of the bosonic field) 

F B 

is an unknown parameter as well. The quantities T", T™, n = 0, 1 (the diagonal components of the 

F B 

energy-momentum tensors), T, T (the corresponding traces), a(ip), A(<p), V(ip), W(a), /(/i), and g(fi) 
are given functions of their arguments, and 7 is a given constant (see below). 

Having in mind the physical assumptions, we have to set the following boundary conditions at 
r = 0: 



A(0) 



dip 
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(0) = 0, a(0) 



fx(o) = He, K R s) = 0, 



(2) 



where a c and fi c are the values of density of, respectively, the bosonic and fermionic matter at the 
star's centre r — 0. The first three conditions in (Q) guarantee the nonsingularity of the metrics and 
the functions A(r), ip(r), a(r) at the star's centre. 

In the external domain r > R s there is no fermionic matter, i.e., one can formally suppose that 

F 

the corresponding density fj,(r) = if x > R s . The fermionic part T™ of the energy- momentum tensor 
vanishes also identically and, thus, the last equation (with respect to the function fj,(r)) in (1) can be 
removed from the system. 

As it is required by the asymptotic flatness of space-time, the boundary conditions at the infinity 

are: 



v(oo) — p(oo) — a (00) — where 



(•)(oo) lim (.)(r) 



(3) 



We seek for a solution [A(r), v{r\, ip(r) , air) , /i(r) , R s , O] subjected to the nonlinear system (1), satis- 
fying the boundary conditions (0) and (Q). At that, we assume the function fj,(r) is continuous in the 
interval [0, R s ], while the functions A(r), v(r) are continuous and the functions p(r), a(r) are smooth 
in the whole interval [0, 00), including the unknown internal boundary r = R s . 

The formulated boundary value problem (BVP) is a two-parametric eigenvalue problem with respect 
to the quantities R s and f2. 

Let us emphasize that a number of methods for solving the free-boundary problems are considered 
111 detail in |g §. 

Here we aim at applying the new solving method to the above formulated problem. This method 
differs from the one proposed in Q and for the governing field equations written in the forms (1) it 
possesses certain advantages. 
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2 Method of solution 



At first we scale the variable r using the Landau transformation x = r/R s , x £ [0, 1] U (1, oo) ||]. 

For given values of the parameters R s and f2 the independent solving of the system (1) in the internal 
domain x £ [0, 1] (inside the star) requires seven boundary conditions. At the same time, we have at 
disposal only six conditions of the kind (||). In order to complete the problem, we set additionally one 
more parametric condition (the value of one from among the functions \i(x), fi(x), (fi(x), <p'i{x), <Tj(x), 
or <j'i(xj) at the point x = 1). Let us set for example 

= ip Sl ip a — parameter. (4) 

Obviously, the solution y; = {A.;, i^, yjj, <jj, ji} in the internal domain x G [0, 1] depends not only on the 
variable x, but it is a function of the three parameters R s , O, (f s as well. 

Concerning the solution y e = {A e , v e , tp e , a e } in the external domain x > 1 six boundary conditions 
are indispensable for the solution of the equations (1). At the same time, only the three boundary con- 
ditions (|^) are given. Let us consider that the solution in the internal domain x G [0, 1] is knowledged. 
Then we postulate the rest three deficient conditions to be the continuity conditions at the point x = 1. 
The first of them is similar to condition (^) and the other two we assign to two arbitrary functions 
from among \{x), v{x), f'(x), cr(x), and cr'(x), for example A e (l) = A;(l); </? e (l) = fs\ 0e(l) = cr«(l)- 

Let us suppose that the internal and external solutions are known. Generally speaking, for given 
arbitrary values of the parameters R s , f2, and ip s the continuity conditions with respect to the functions 
v{x), (p'(x), and a'(x) at the point x = 1 are not satisfied. We choose the parameters i? s , fi, and ip a in 
such manner that the continuity conditions for the mentioned functions to be held, 

i.e., v e (l, R s , O, (f s ) - 1^(1, R s , fi, (p s ) =0, 

<p' e (l, R s , n, Vs ) - ^(1, R„Sl, <p.) =0, (5) 

a'ji, R s ,n, <p a ) - ^(i, R s ,n, <p a ) =o. 

These conditions should be interpreted as three nonlinear algebraic equations in regard to the unknown 
quantities R s , Q and ip a . 

So, we have to solve the nonlinear BVP (1) - (||) in conjunction with the algebraic system (5). The 
traditional technology for solving such problems, based on the methods like shooting ||, has some well 
known disadvantages. 

In this paper we use the Continuous Analogue of Newton Method for solving the above nonlinear 
spectral problem which proposes a common treatment of both differential and algebraic problems. 
Detailed description of the algorithm can be found in jj] . 



3 Some numerical results 



For the purpose of illustrating we will shortly discuss some results obtained from numerical experiments. 
In the present article we consider concrete scalar-tensor model with functions (see Section |l|) 



A(ip) =exp 



i 



V{v) = -[1-A\ v )]\ W(a) 



/(/i) = - [(2fx - 3) V/i + + 3 In (yjl + ^TTJl 
<?(M) = \ [(6// + 3)V/i + /i 2 -31n y/1 + /i) 



The quantities b and A are given parameters. The functions f([i) and g(/x) represent in parametric 
form the equation of state of the non-interacting neutron gas, while the function Wicr) describes the 
boson field with quadratic self-interaction. 

The behaviour of the calculated eigenfunctions cr(x), tp(x), v{x), and fi(x) is typical for wide range 
of the parameters and was discussed in detail in |Q. 

From a physical point of view, it is important to know the mass of the boson-fcrmion star and the 
total number of particles (bosons and fermions) making up the star. 
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Parameters: A = 0, y = 0.1, b = 1, o = 0.002 
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Figure 1: The star mass M and the rest fermion mass Mrf as functions of the central value fi c . 



The dimensionless star mass can be calculated by the formula: 



M= r 



,.2 



dr. 



The dimensionless rest mass of the bosons (total number of bosons times the boson mass) and the 
dimensionless rest mass of the fermions are given by: 

oo oo 

Mrb = SI J r 2 A 2 {tp) exp ° 2dr , M RF = b J r 2 A 3 {p) exp (j-) n(/x)dr, 

o o 

where n(ii) is the density of the fermions. In the case we have = /ji(x). 

The dependencies of the star mass M and the rest mass of fermions Mrf on the central value 
fi c of the function (i(x) are shown in the configuration diagram on Fig. 1 for A = 0, 7 = 0.1, 6 = 1 
and er c = 0.002. It should be pointed that for such small central value a c we actually have a "pure" 
fermionic star. On the figure, it is seen that from small enough values of /i c , in the vicinity of the the 
peak the rest mass is greater than the total mass of the star, which means that the star is potentially 
stable. 

On Fig. 2 the binding energy of the star E\, = M — Mrb — Mrf is drawn against the rest mass of 
fermions Mrf for A = 0, 7 — 0.1, 6 — 1 and a c = 0.002. Fig. 2 is actually a bifurcation diagram. By 
increasing the central value of the function fi{x) a cusp appears. The presence of the cusp shows that 
the stability of the star changes - one perturbation mode develops instability. Beyond the cusp, the 
star is unstable and may collapse, eventually forming a black hole. The corresponding physical results 
for pure boson stars are considered in our recent paper ||. 
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